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We propose a biochemical model providing the kinetic and energetic descriptions of the processivity dynamics 
of kinesin and dinein molecular motors. Our approach is a modified version of a well known model describing 
kinesin dynamics and considers the presence of a competitive inhibition reaction by ADP. We first first recon- 
struct a continuous free-energy landscape of the cycle catalyst process that allows us to calculate the number 
of steps given by a single molecular motor. Then, we calculate an analytical expression associated to the 
translational velocity and the stopping time of the molecular motor in terms of time and ATP concentration. 
An energetic interpretation of motor processivity is discussed in quantitative form by using experimental 
data. We also predict a time duration of collective processes that agrees with experimental reports. 
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I. INTRODUCTION 

Molecular motors are ubiquitous entities in living cells 
that participate in controlling and synchronizing many 
cellular transport processes and are involved in different 
metabolic pathways . This essential role makes the study 
of their biochemical kinetics of relevant interest in biol- 
ogy and medicine 2 . In addition, the cyclic operation of 
these motors and its possible control may be of help for 
developing future externally controllable nano-machines 

3,4 

Many experimental and theoretical studies have been 
devoted to understand the detailed biochemical reactions 
associated to the motion of kinesins and other molecular 
motors, see, for instance, Refs. 3 ' 5-12 . One of the aims 
of these studies is to characterize protein conformational 
changes in order to understand the sequence of reactions 
accompanying motor translation and determining its en- 
ergetic dependence 6 " 8 ' 13 . To this end, single-molecule 
fluorescence spectroscopy and polarization microscopy 
techniques have been used in order to better discriminate 
individual motor activities due to the large heterogene- 
ity in motor behavior . These studies allow an accurate 
determination of step sizes, rates and pausing properties 
of different motor dynamics that can be used to propose 
more specific theoretical models 6 ' 9-12 ' 14-16 . Often, these 
models follow two different approaches, one describing 
the collective behavior and properties by viewing the mo- 
tors as catalytic agents of a chemical reaction 6 ' 14 whereas 
the second one focuses on the individual description that 
analyzes the specific conformational changes occurring 
during motility 9-12 . From a more physical point of view, 
several models taking into account the stochastic and 
thermodynamic nature of the translational dynamics of 
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molecular motors were proposed and used to describe 
some general features of molecular motor activity 15-21 . 
The essential ingredient to perform the description in 
these approaches is to provide a model for the energy 
landscape of the cyclic operation of motors 15-17 ' 20 ' 21 , 
however these models are barely related to the chemical 
reactions that drive the process. 

Despite of the particular mechanisms determining the 
rate at which each catalyst cycle occur, some general 
characteristics can be used in order to reconstruct the 
energy landscape of the entire process, that is, they allow 
for determining the processivity of the molecular motor 
under consideration. Understanding the kinetic and en- 
ergetic properties associated to the processivity of molec- 
ular motors is of central importance when connecting the 
specific dynamics to the biological function ' . This is 
because cellular transport processes are usually driven by 
more than one motor, therefore requiring a high degree 
of cooperativity. 

As proteins, enzymes and molecular complexes, the 
processes developed by molecular motors use chemical 
energy stored into molecules such as ATP or GTP, pro- 
duced by the mitochondrial system of the cell 2 . However, 
many studies are performed in vitro with well controlled 
ATP concentrations and somehow simplified conditions 
concerning the viscoelastic properties of the surround- 
ings. This allows for a better analysis of the detailed dy- 
namics of the motor but may hide some aspects of its per- 
formance in vivo, such as the transport of proteins, RNA, 
vesicles and even organelles 23 that may be related to, for 
instance, several exocytosis-endocytosis processes 24-27 . 

In this work, we analyze the thermodynamics and 
chemical kinetics of two theoretical models describing 
motor processivity 6 ' 9 ' 28 , their time activity and the asso- 
ciated energy consumption. These two hand-over-hand 
models were first developed to describe the particular 
problem of intracellular transport via kinesin along mi- 
crotubules and are based on the well known evidence that 
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the energy used for the motion comes from ATP hydroly- 
sis 6-9 . We first simplify the general schemes proposed in 
the literature by assuming that the three initial steps of 
the reaction sequence are slower than the last ones, and 
therefore determine the time step of the overall reaction. 
In a second model, we also assume that the ADP pro- 
duced during the motion of the molecular motor plays 
the role of a competitive inhibitor 9-12 . In both cases, we 
use the concept of the degree of reaction 29 and follow the 
rules of thermodynamics to reconstruct the correspond- 
ing Gibbs free-energy landscape for a single cycle and 
motor. In this way, the free-energy we obtain is a quan- 
titative representation of an enzyme catalytic reaction 
in which protein conformation fluctuations in the pres- 
ence of the substrate are taken into account, that is, it 
constitutes a model of an enzyme complementary to the 
transition state and not to the substrate . 

Then, using Fourier analysis we are able to recon- 
struct the complete energy landscape of the whole cat- 
alytic reaction and the translational velocity of the mo- 
tors through a well established procedure for the analysis 
of the enzymatic reaction scheme. By its nature, this sec- 
ond information is important because complements that 
of the free-energy by giving a collective notion associated 
to motors activity. In the case of the process with inhibi- 
tion, we also determine the time course of both, the free 
energy and the translational velocity. We analyze the ef- 
fect of different ATP and ADP levels on this velocity. In 
addition, we also give a simple expression for the stopping 
time of a collectivity of molecular motors, an important 
quantity characterizing the finite processivity of the mo- 
tors when participating in metabolic tasks. These results 
allow us to predict the number of steps that a single mo- 
tor may perform before stopping, and also the average 
traveled distances and the associated ATP consumptions 
of a collectivity of motors giving the initial values of the 
different parameters entering in the description. The ef- 
fect of the load is also taken into account. 

The article is organized in five sections. In Section II 
we present the essential biochemical reaction model that 
we will modify accordingly with the aim to determine, 
in Section 3, the continuous free-energy landscape of the 
catalytic reaction and the average translational velocity 
of a kinesin motor. Section IV is then devoted to ana- 
lyze the inhibitory effect of ADP as the main energetic 
mechanism determining the finite time of motor proces- 
sivity. The time dependence of the translational velocity 
and the stopping time are also calculated and discussed 
using values for the experimental parameters taken from 
literature. 



II. THE HAND-OVER-HAND MECHANISM 

The two models accounting for the transport pro- 
cess via kinesin activity that we will discuss in the fol- 
lowing sections are simplifications based on the orig- 
inal biochemical reaction scheme proposed in Ref. 



and based on detailed experimental studies citeViss- 
cherl999,ADP,svoboda,hanckock. The first reaction of 
a sequence of six (as shown schematically in Figure 1) 
accounts for the capture of an ATP molecule (T) by the 
microtubule-kinesin complex MK in order to produce 
the enzyme-substrate complex MKT a 

MK a +T ^ MKT a , (1) 

where the subindex a indicates that the corresponding 
kinesin head K a is attached to the microtubule, whereas 
the second head /3 is free. According to this scheme, the 
complex MK plays the role of an enzyme that acts over 
substrate T through a catalytic reaction. This means 
that free kinesins cannot consume ATP periodically in 
time unless they are attached to a microtubule . The 
complete reaction sequence associated to a single step is 
given through Eqs. (2)- (6) 

MKT a ^ MK^ , (2) 



MKT' a ^ MKT a Dp , (3) 



MKT a Dp ^ MKpT a , (4) 



MKpT a — ^> MKp(D o P) a , (5) 



MK P (D o P) a % MKfs + D + Pi . (6) 

In accordance with experimental observations 6 ' 8 ' 28 , the 
last reaction reflects the fact that the hydrolysis of ATP 
at the active site a of the kinesin produces enough en- 
ergy in order to liberate the corresponding head. After 
this reaction occurs, the cycle is completed and an initial 
state MKp is recovered one step forward from the initial 
position and with lower free-energy ref. . The number 
of repetitions of this cycle that the molecular motor is 
capable to perform determines its processivity. 

Here, it is important to stress that reaction (6) pro- 
duces an ADP molecule and an inorganic phosphate Pi 
that can also react with the corresponding active head 
of the kinesin . As a consequence of this, a complex 
inhibition scheme controlling motor processivity may 
arise 9 ' 30-32 . Since during the evolution in time of the 
process ADP and Pi concentrations grow, then the proba- 
bility of occurrence of an inhibition event increases. This 
process and its influence on motor processivity will be 
analyzed in detail in the following sections. 
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(3) 

MKT' a i MKT a D p 



(4) 

MKT a D B % MKrTu 








H ■ 




■ □ 




(6) 



o P) a % MKp + D+ 




FIG. 1. Schematic representation of Eqs. (l)-(6). Blue molecules represents ATP, yellow molecules represents ADP, purple 
molecules inorganic phosphate. 



III. GIBBS FREE-ENERGY LANDSCAPE 
RECONSTRUCTION 



In order to reconstruct the free-energy landscape asso- 
ciated to the motion of the molecular motor as modeled 
by Eqs. (l)-(6), we will first assume that reactions (4)- (5) 
are much faster than reactions (l)-(3) 6 ' 7 , and therefore 
its contribution to the total time of the global reaction 
can be neglected. As previously mentioned, reactions (1)- 
(3) are associated to the formation of the active enzyme 
MK by fixation of kinesin a-head to the microtubule 
in the presence of an ATP molecule, reaction (1). This 
process produces the formation of the enzyme-substrate 
complex MKT a that may fluctuate allowing a conforma- 
tional change of the dimer constituting the motor stalk 
and producing the secondary enzyme-substrate complex 
MKT' a1 reaction (2). The secondary complex has the 
appropriate structural conformation in order to attach 
the /3-head to the microtubule and rapidly hydrolyze the 
ATP of site a in order to promote the translational mo- 
tion, reactions (3)- (6). According to this model, the iso- 
merization reaction (2) is the key ingredient in order to 
reconstruct the Gibbs free energy landscape and implies 
the presence of a small energy barrier that can be ther- 
mally overcome. 

If we assume that the reactions taking place during 
the motion of the molecular motor are not far from equi- 
librium, we may use non-equilibrium thermodynamics 
to analyze the chemical kinetics of the process 29 . For 
each elementary reaction, the differential variation of the 
Gibbs free-energy in terms of the degree of reaction £ can 



be written in the form 

dGj = — Uijfiijd^j , 



(7) 



ZJ and Hij are the stoichiometric coefficient and 



where 

the chemical potential of the i-th species in the j-th reac 
tion. The relation between the mass fraction rii and the 
corresponding degree of reaction £j is: drti = —Vijd^j, 29 



Thus, assuming that the system is sufficiently diluted we 
have 



dGj = -RT^Vij In 



(8) 



where we used the expression: Hi = RT]n\rii/n^ q \, with 
rti the molar fraction of the z-th chemical specie and n^ q 
its equilibrium value. 

In this form the change of Gibbs free-energy of reaction 
(1) can be written in the form 



dG 



-RT In 



MK 



MKT 



MK e 



^MKT 

MKT et 



df. (9) 



Here, the relation between £ and the corresponding molar 
fractions is established in Table I for a time r elapsed 
after the chemical reaction started. 

Now, taking into account reaction (2), Eq. (9) can be 
rewritten by using the fact that the complexes MKT and 
MKT' are equilibrated according to reaction (2), see the 
Appendix A for details. As a consequence of this, they 
obey the relation: = umkt' I^mkt-, where K^ = 
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Time 



Molar fractions 



t = n° MK n°MKT n °MKT' 

t T n° MK - £ n° T - £ n° MKT + £ ^°mkt' ~ € 

TABLE I. Stoichiometric ratio of chemical species. From this 
table it follows that for reaction (1): vmk — 1, vt — 1 and 
vmkt = — 1. The last column takes into account the fact 
that the equilibrium reaction (2) also produces MKT, and 
therefore v MKT i — 1. The superscript ° indicates the initial 
values of the mass fraction the corresponding species. 

K Q e~^ ] / kBT is the equilibrium constant corrected by the 
factor e -f 6 / k sT which takes into account the effect of the 
load / (weight) imposed by the cargo on the motor and S 
corresponds to the distance travelled by the cargo during 
the isomerization , typically one half of the characteristic 
step distance of 8 nm. Thus, for cargo dependent motions 
we have in turn 



eq 

MKT' 



MKT 



K n e~f 5 / k B T 



(10) 



Substituting now the previous relation into (9), using 
Eq. (A8) obtained in the appendix and performing a 
Taylor expansion of the logarithm: ln(x) ~ x — 1 for 
x ~ 1 (notice that the condition x = rii/n^ q ~ 1 is 
always obeyed for reactions not far from equilibrium) we 
finally obtain the following expression for the effective 
Gibbs free-energy change of reactions (1) and (2) 



W = - [(n°MK ~ OK " 0(n° MKT , - 



(11) 



-(« 



eq eq 
MK T 



i% KT ,)+f5/k B T}dt. 



An integration over £ yields the following expression 
for the Gibbs free-energy change AG(£) 



AG(C) 
RT 



=Ae+Be+ce+{D-D eq )t 



+f5/k B T£, 



(12) 



where we have defined the constants 

I u — n °MK + n ' 

4 

1 



A 



B 



l MKT' 



C — 2 ( Tl MK n T + n °MK n °MKT 
D = ^MK n T n °MKT' 5 D eq 



(13) 



°MK'°T ,l MKT> • 



It is important to emphasize here that coefficients B, C 
and D depend upon the initial concentrations of enzyme 
MK. ATP and MKT whereas D eq depends upon the 
equilibrium concentration values of these species. 

Equations (12) and (13) constitute the basis for re- 
constructing the free-energy landscape of the catalytic 
reaction and therefore are useful to explain the proces- 
sivity of a single molecular motor. Figure 2 shows the 
normalized Gibbs free energy (12) in terms of the degree 




FIG. 2. The normalized free energy difference AG in terms of 
the reaction coordinate £. The blue solid line corresponds to 
the case without load (cargo weight) whereas the green dashed 
line corresponds to the case when a cargo of / = 0.02pN is 
applied to the molecular motor. The differences of the free 
energy between states depend on the initial concentrations 
of enzyme MK, ATP and MKT' whereas the global tilt- 
ing of the potential depends on the equilibrium conditions 
incorporated by D eq . The cargo weight reduces the free en- 
ergy difference between the isomeric states MKT and MKT 1 , 
thus decreasing the average translational velocity of the mo- 
tor. The initial values used are: n MKT , = 0.02, n° T = 0.78 
and n° MK = 0.2. The values of the rate constants are given in 
Fig. 6. 



of reaction £ for values of the parameters taken from ex- 
periments . The free-energy difference is tilted and 
is asymmetric because the left maximum at (£ = 0) is 
higher than the right one (£ = 1). In addition, the MKT' 
state is energetically more favorable than the MKT one, 
fact that promotes the motion of the motor. In the fol- 
lowing sections, we will show that the present scheme for 
the catalytic reaction can be generalized by considering 
competitive inhibition by the final products of the se- 
quence of reactions (the production of ADP and Pi) and 
how this leads to invert the above situation, that is, when 
due to the presence of inhibitors the MKT state becomes 
energetically more favorable than the MKT' one, which 
results in the end of the catalytic reaction. In physi- 
cal terms, these asymmetries shown by the potential im- 
ply that the catalytic free-energy landscape we deduced 
from the biochemical reaction scheme is a type of poten- 
tial called "tilted ratchet" potential 15 " 17 ' 19 ' 21 . Thus, the 
biochemical model proposed here gives an experimental 
basis to these type of potentials heuristically proposed in 
the literature. 

In Figure 2, the local maximum at £ ~ 0.54 corre- 
sponds to the free energy barrier associated to the iso- 
merization reaction (2), also represented in the schematic 
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view of the motor on the top of the figure. The corre- 
sponding activation energy will be denoted by and the 
difference between isomeric states MKT and MKT' will 
be denoted by AG^ (see also Figure 2). The solid blue 
line corresponds to the motion of the kinesin with a cargo 
of negligible weight. The green dashed line shows that 
the tilting of the potential decreases when the weight of 
the cargo is increased. The values of the parameters used 
to represent AG were taken from Refs. 7 ' 9 . 

A. Processivity: The single motor free-energy landscape 

The free-energy difference (12) only applies for describ- 
ing the isomerization in a single step, and not the over- 
all process associated to the displacement of the motor. 
Thus, in order to model the whole catalytic process, first 
we have to consider that it constitutes a sequence of sin- 
gle steps, one after one. 



T 




T T T 



1 2 3 t 



ATP — > ADP + P % 

reaction coordinate 



that the reaction occurs with constant initial values of 
the mass fractions n° MK , n° T and n° MKT , as illustrated in 
Table I. 




FIG. 4. First three steps of the free-energy landscape associ- 
ated to kinesin translation as obtained from Eq. (28) for the 
initial values: n° MKTl — 0.02, — 0.78 and n° MK — 0.2. 
Two energy barriers appear, one related to the isomerization 
reaction (2), G\ and other associated to the a — /3-shifting 
for complex MKT' a to MKT/3, G a -p which is surmounted via 
ATP hydrolysis. The inset shows the free-energy differences 
AC7 f between the states MKT a and MKT' a , and between the 
sates MKT' a and MKTp, AG a -(3- The values of the other 
constants were taken from' and and are given in Fig. 6. 



FIG. 3. Graphical interpretation of the relation between the 
degree of advance of a single reaction £, the time duration 
of each reaction r and the global degree of reaction £ g corre- 
sponding to the catalytic reaction. Each step is characterized 
by the hydrolysis of one molecule of ATP. 



Following this idea, the catalyst landscape (cyclic-steps 
landscape) can be viewed as a periodic potential whose 
elemental cycle is described by Eq. (12). In order to 
perform the catalyst landscape reconstruction, we shall 
generate a periodic potential starting from Eq. (12) and 
performing a Fourier series expansion considering only 
the terms depending on the initial values of the molar 
fractions A, B, C and D, see Eq. (13). The other two 
terms related to D eq and /5/KbT are exclusively deter- 
mined by the equilibrium conditions of the problem and 
of the load / applied by the cargo. 

Hence, the catalyst energy landscape can be expressed 
in terms of the global degree of reaction £ g that goes from 
to 00 (as the time variable r) and which is related with 
the reaction coordinate £ as it is shown in Figure 3. Each 
cycle (step) has a duration of r ~ 10 ms and the first step 
goes from to r whereas the second step goes from r to 
2t, etcetera. In the present case, the lack of inhibition 
by the product of the catalytic reaction allows to assume 



Performing the Fourier series of Eq. (12), the resulting 
expression for the catalyst energy landscape in terms of 
the global reaction coordinate is 



1 

RT 



cyclic 



A 



■ ^ B n sin(27rfc£ g ) - D eq £ g 



k B T 



(14) 



Co, 



in which the Fourier coefficients A n and B n are defined 

by 



A k = [ {A£ 4 + + Ci 2 + D£) cos(27rfc£)d£, 
Jo 

B k = I { Ai A + Bi z + d 2 + DO sin(27rfc£K , 
Jo 



(15) 



In order to obtain Eq. (14) we have assumed that, for 
the first cycle, the equilibrium values of are equal to 
the initial values n°, thus implying that D = D eq . The 
first three steps of a kinesin motor are represented in 
Figure 4 via the Gibbs free-energy landscape (14) of the 
catalytic reaction. An interesting feature of the model is 
that appears a second energy barrier G a -$ (see Figure 
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4) that separates the state MKT^ from state MKTp [see 
also Eqs. (3)-(6)]. The amplitude of this barrier is ten 
times larger than that corresponding to the isomeriza- 
tion reaction. The fact that this second energy barrier is 
not too high, seems to agree with the observation 26 that 
some kinesin motors may operate in two forms, one with 
directed motion by using ATP and a second non-directed 
motion without using ATP but using the thermal energy 
at their disposal. Finally, the Gibbs free-energy differ- 
ence between the states MKT' a and MKTp is denoted 
by AG a - p. Numerical values for the free-energy bar- 
riers and differences were estimated using experimental 
data 7 ' 9 for the initial concentrations in Eq. (14) and are 
presented in Table II. Since the initial values of n° MK , 
n9p and Ti° MKT , are constant, no change of the amplitude 
of the barriers and energy differences appears during the 
catalytic process. This clearly represents an idealization 
that will be improved in following sections were the effect 
of ADP inhibition will be taken into account. 



Energy kJ/mol Step 1 Step 23 Step 62 

G T 0.012 0.012 0.012 

AG f -0.017 -0.017 -0.017 

G a -/3 0.083 0.083 0.083 

AGa-p -0.005 -0.005 -0.005 

TABLE II. Magnitudes of the energy barriers and differences 
associated to the free-energy landscape of the catalytic reac- 
tion as obtained from the model (14) using the same values 
as those used in Figure 4. 



B. Motor translational velocity 



Menten approach is sufficiently accurate for large con- 
centration values of ATP ' . 



IV. QUANTITATIVE INFLUENCE OF ADP ON 
PROCESSIVITY 

Experimental studies have reported that the presence 
of ADP and Pi reduces the processivity of kinesins by 
inhibiting the formation of the enzyme-substrate com- 
plexes MKT and MKT' ', which are essential in turn to 
the motility of the motor 9 ' 30-32 . On the basis of these 
observations, the previous analysis can be modified ap- 
propriately in order to account for this finite processivity 
of kinesins, that is, their limited number of cycles. 

Several types of ADP and Pi inhibitions are possible 
in principle, like competitive or non-competitive depend- 
ing on their release order in the corresponding kinesin 
head . Each type of inhibition will modify in different 
way the expression for the kinesin velocity. The explicit 
form in which these corrections will appear are described 
in detail in Appendix D. However, it was experimentally 
found in Ref. 9 that the effect of the non-competitive and 
competitive inhibitions by Pi are only relevant at low 
ATP concentrations, which suggests that they could be 
observable only for a short time period at the end of the 
process. This means that competitive inhibition by ADP 
can be considered as the leading inhibition effect of the 
process, as it will be considered here. The generaliza- 
tion of the model to include the other type of inhibition 
mechanisms is straightforward. 

In order to deduce the free-energy landscape contain- 
ing the effect of ADP inhibition, it is necessary to first 
analyze how motor displacement velocity changes, and 
how the concentrations of ATP and ADP evolve in time. 



It is well know that kinetics of most part of biologi- 
cal processes mediated by non-allosteric enzymes are well 
described by Michaelis- Menten like equations . This fact 
may be used in the mechanism scheme (l)-(6) by as- 
suming that reaction (3) is the slow step in reaction- 
mechanism, as we did in the previous subsection. Thus, 
assuming that steps (4)- (6) are fast , the translational 
velocity v([T]) can be obtained by multiplying the reac- 
tion velocity w([T]) by kinesin's step distance, d = 8nm, 
and divided by the total enzyme concentration [MK] Q1 
which measures the number of active motors in a given 
time. The final expression is 6 



v([T}) 



dk cat [T] 



K 



[T] 



(16) 



where k cat = k^K^ l{K^ + 1) is the catalyst rate con- 
stant, denoting the maximum number of enzymatic reac- 
tions catalyzed per second 6 . K^ M is the pseudo-Michaelis- 



Menten constant defined by K\ 



A more detailed derivation of Eq. (16) is reported in 6 . 
Experimental results have shown that the Michaelis- 



A. Motor translational velocity in the presence of ADP 
inhibition 

In accordance with the previous considerations, we will 
assume that ADP inhibition can be described as a com- 
petitive reaction since kinesin active sites seem to be 
affine to ATP or ADP in an exclusive way. This hypoth- 
esis modifies the reaction scheme (l)-(6) by including the 
parallel reaction 



MK 



D ^ MKD„ 



(1-1) 



Following the same procedure as in section 3.2, it may 
be shown that the displacement velocity of the kinesin in 
the presence of inhibition is 



v([T],[D}) = 



dk cat [T] 



(20) 



where Kj is the inhibition constant . It is important to 
emphasize that this relation assumes step (3) as the slow 
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MK a + D MKD 







a u 




■ 


□ 




FIG. 5. Schematic representation of parallel inhibition reac- 
tion, Eq. (1-1). 



step that determines the behavior of the kinetics, but it 
does not takes into account that ADP is a product of 
the whole reaction scheme. As a consequence of this, Eq. 
(20) can not be used to evaluate the ADP production 
kinetics, this is more probably determined by the last 
step (6). 

In order to reconstruct the free energy landscape of the 
catalytic process in presence of inhibition by ADP, it is 
necessary to implement the time dependence of ATP and 
ADP concentrations since in the overall process the initial 
concentrations after each cycle (n^, n Q MK and n Q MKT ) 
change as a function of time. 

For this purpose, it is necessary to solve the reaction 
velocity equation w — —d[T]/dt. The solution is given in 
implicit form by the expression 



K ] In 



[T]o 



Pi* - mo = 



(21) 



where the subindex t indicates time dependence, and 
v max = k cat [MK] . From these relation it follows that 
[T]t is given by the W Lambert function 33 . 

For the values of the parameters previously used, it 
may be shown by directly evaluating the function (or by 
numerically solving the equation) that Eq. (21) can be 
well approximated by 



[T]t^[T] 



(22) 



Here, it is convenient to mention that due to the 
crowded nature of the intracellular medium, transport 
of ATP and ADP by diffusion becomes a very slow and 
probably even a confined process, as it has been observed 
experimentally for chromosomes and theoretically ex- 
plained in . As a consequence of this the theoretical 
treatment of the process can be performed, in a first ap- 
proximation, without taking into account spatial hetero- 
geneity or ATP production. 

On the other hand, according to experimental 
results 10 ' 11 and the previous considerations on Eq. (6), it 
follows that the time behavior of the ADP concentration 
can be modeled by the first order kinetics of the form 




FIG. 6. Translational velocity of a kinesin motor as a func- 
tion of time given by Eq. (24) with the following val- 
ues for the constants entering into Eqs. (24) and (25): 



[T] ImM, [D] 0, v x max 9UU nms \ v ma x 
11.25/xMs" 1 , [MK] lOOnM, K ] M 28/xM, k cat 
113s -1 , fei = 4.0/iM- 1 s-\ k-i = 3.8 x 10" 5 s"\ k 5 = 
5.6 x 10 3 s~\ K = 0.0017 were taken also from Refs. 7 ' 9 . 
The value K X M = 2.24 fim was calculated using its definition 
Kl = dKl/[MK] . 



Using equations (22) and (23) we obtain the following 
expression for the average displacement velocity in terms 
of time 



v(t) = 



V max(ts t) 

h + (t s - 1) ' 



(24) 



where we have defined the maximum displacement ve- 
locity = dk cat = dv max /[MK] , the stopping time 
t s = \T] /v max and the function 



Vrnax 



(25) 



where K^ M 



[D] t = [D} ^[T} (l-e-^). 



(23) 



dK ] M /[MK] and [D] t is given in (23). Eq. 
(24) is shown in Fig. 6 for different values of the inhi- 
bition constant Kj. It is notable that when inhibition 
is small (Kj large) the displacement velocity is almost 
constant during the elapsed time of the process and near 
to the saturation velocity. When inhibition increases (Kj 
decreases) the velocity departs from saturation levels and 
becomes non constant. These result may be important 
for biological processes such as exocytosis 36 . In Fig. 6 
we have assumed that the process starts instantaneously 
when there exist saturation levels of ATP and Ca + 2 sig- 
nalment have occurred . We have modeled this situation 
by means of a Heaviside function at a initial time U in 
Eq. (24): v(t) = Q(t-U)v(t). 

We can estimate the time duration of a given process 
using Eq. (24) with the values used on Fig. 6. We obtain 
t s = [T] /v max = 400s. Because the stopping time has 
been calculated by assuming that the translational ve- 
locity depends on the total concentration of the enzyme 
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MK, it contains the effect of motor cooperativity and 
has to be compared with cellular processes that prob- 
ably involve the participation of several motors. This 
may be the case of some exocytosis-endocytosis process. 
In fact, the stopping time we obtained agrees well with 
the characteristic times observed in experiments study- 
ing membrane resealing, where the active transport of 
vesicles from the inner part of the cell to the plasma 
membrane is very probable (see Ref. ). Another impor- 
tant result emerging from the present analysis is related 
to the fact that the constant K^ M = cIKm/[MK]o can 
be interpreted as the average distance of advance un- 
til its velocity down to half of saturation velocity ^ aa .. 
The average travelled distance corresponding to the es- 
timated stopping time is K} M ~ 2.4/xra, thus implying 
the consumption of about ~ 300 ATP molecules per mo- 
tor. These results also compare well with experimental 
evidence 28 . 



B. Single motor free-energy landscape in the presence of 
ADP inhibition 

The analysis and results obtained in previous sections 
bring us key facts to incorporate ATP consumption and 
ADP production into the free-energy model. As we 
showed, the increase of ADP concentration may modify 
the equilibrium on (2) because a parallel inhibition reac- 
tion between the MK complex and ADP is also present 9 . 
However, one may still assume that free enzyme is in a 
steady state in similar way to the Michaelis-Menten ki- 
netic 6 ' 9 ' 12 ' 37 ' 38 . 



Time 



Molar fractions 



n° T (t) 



'(*) 



t + r n° MK -i n° T {t + r)-^ n° MKT , (t + r) - £ 

TABLE III. Time dependent stoichiometric ratios. The molar 
fraction of MK is independent of time since it is considered 
the free enzyme. 



The procedure to obtain the single motor catalytic en- 
ergy landscape in the presence of ADP as an inhibitor is 
similar to the one followed in Section 3. However, in the 
present case we have to consider that after each single cy- 
cle, the initial values for the mass fractions n° MK , n?p and 
n °MKT' nave to be recalculated by taking the correspond- 
ing values n^Jr), u t( t ) anc ^ u mkt / ( t )- Thus, for any 
other subsequent cycle, the data of Table I become: n°(t) 
and n?(t + r) - f , see Table III. 

This time dependence enters because the initial con- 
centration of ATP decreases with time r as the catalytic 
reaction takes place and produces ADP, Eqs. (22) and 
(23). Hence, the coefficients A n , B n in (14) become func- 
tions of time r through the initial concentrations after 



each cycle, n°(r). This may be represented schemati- 
cally by the subindex r in coefficients: B r , C T and D T 
(see Appendix B for details). 

The explicit form of this time dependence can be de- 
rived following Ref. . It may be shown (see Appendix 
B) that ATP concentration decreases following the re- 
lation: n^i^g) = n9p — kcatr^g, which can be obtained 
from Eq. (22) and considering the existing bijection be- 
tween r and £ g indicated in Figure 3. In addition to this, 
the parallel inhibition reaction between the MK com- 
plex and ADP (ADP competes for the same active site 
of kinesin as ATP) avoids the formation of both enzyme- 
substrate complexes MKT and MKT. Thus, this inhi- 
bition increases with £ g because the whole reaction cycle 
produces ADP and therefore reduces MKT'. This de- 
crease of MKT' cannot be directly inferred from the ki- 
netic equation since the reaction velocity do not depends 
on the mass fraction of the complex MKT' . However, in 
similar way as estimated for the ATP decrease, we can 
assume that the initial mass fraction rt° MKT , decreases in 
linear form with that is: n° MKTI {^ g ) = n° MKT , -K Q £ g , 
where the constant K Q is the equilibrium constant with- 
out load [see Eq. (2)]. 

Thus, after performing the Fourier series having in 
mind the previous considerations, the resulting expres- 
sion for the catalyst energy landscape is 



1 A (t) 

AG^ g: r) cycHc ~ + J2 A k (T)cQB(2irk£ g ) 



RT 



2=1 



2=1 



- B k (r) sin(27rfcy - D eq ^ g + , 



(26) 

in which the Fourier time dependent coefficients An(r) 
and B n (r) are defined by 

A k (r)= [ (AZ 4 + B r f + C r Z 2 +D r Z)cos(27rkZ)d^ 
Jo 

Bjb(r) = / (A^ + B T f + C T e + D r £) sm(2nk^ . 
J o 

(27) 

Establishing now the bijection between r and £ g ac- 
cording to Figure 3, we may recast Eq. (26) in its final 
form 



1 n(t \ ~ ^ofe) 

^\S,g)cyclic — 



RT 



J2 Ak(€g)cQa(27rk£g) 



2 = 1 



-^B k (Qsm(2nkQ - D eq £ g 



2=1 



fS 
knT 



(28) 



As in Section 3, in order to obtain Eq. (28) we have 
assumed that, for the first cycle, the equilibrium values 
of rii's are equal to the initial values n?'s, thus implying 
D = D eq . Equation (28) is very general since models 
both, every single step during the translation of a single 
kinesin motor, but also the overall process and even the 
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finite processivity of the motor, as we will explain in the 
next. 



Onm 8nm 16nm 24nm 32nm 40nm 




61 62 63 64 65 66 



Steps 

FIG. 7. Free energy landscape obtained via the Fourier ex- 
pansion Eq. (28) with k = 25 modes for a catalytic reac- 
tion of many cycles with no cargo. The equilibrium values of 
the mass fractions used were: n^ KT , = 0.02, = 0.78 
and = 0.2. a) from cycle 1 to 5 and corresponding nm 

displacement, b) step 33 to 38. c) step 61 to 66. The values 
of the constants used are given in the caption of Figure 6. 



V. DISCUSSION 

Three different sequences of steps of the Gibbs free- 
energy landscape (28) of the catalytic reaction associated 
to kinesin translation are represented in Figures 7 and 8 
for the free of cargo and cargo cases, respectively. In sim- 
ilar way as in Section 3, Eq. (14), in the present case we 
find two energy barriers associated to the translations of 
the kinesin. The first barrier separates the isomeric states 
MKT and MKT 1 , whereas the second one separates the 
a — ^-shifted states. The important fact to emphasize 
here is that the free-energy landscape (28) contains am- 
plitude coefficients A n and B n that depend explicitly on 
the global degree of reaction £ g and the initial mass frac- 
tions n®. Since these dependences take into account the 
feedback inhibition by the product of the reaction ADP, 
which increases with these amplitudes of cosine and 
sinuous terms of the expansion increase after many steps. 
As a consequence of this effect, the magnitude of both 
free energy barriers increases and also changes the tilting 
of the overall potential. 
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Energy kJ/mol Step 1 Step 23 Step 62 

G ?t 0.012 0.018 0.023 

AG 1 " -0.017 -0.014 0.009 

G a -f3 0.083 0.088 0.135 

AGa-p -0.005 -0.07 -0.02 

TABLE IV. Comparison of energy requirements and energy 
differences in three different steps of the catalytic reaction as 
predicted by Eq. (28). These energies were obtained for the 
same values as those used in Fig. 7. 

These two features of (28) have the consequence that, 
after a finite number of steps, the displacement of the 
motor stops. Specifically, this occurs because the fast 
equilibrium between MKT and MKT' states is inverted 
in the sense that the backward state MKT more favor- 
able than the forward state MKT' \ see the last steps in 
Figures 7c and 8c. In accordance with the energy min- 
imizing statement, this change of relative values of the 
corresponding free-energies implies that the motors can- 
not follow in its motion. In addition, the amplitude of 
both barriers increases as the reaction takes place. This 
fact slows down the activated process between MKT and 
MKT' sates and finally makes more improbable that the 
energy coming from ATP hydrolysis is enough to over- 
come the barrier of the a — /3-shifted states, see the dif- 
ference between Figures 7a and 7c, and of Figures 8a 
and 8c. Using the data of the previous analysis, we pre- 
dict that the free of cargo motion is characterized by a 
processivity of about 66 steps, which corresponds to an 
energy consumption of 66 ATP molecules and a traveled 
distance of 528 nm. This number suggests that a mo- 
tor traveling 2.4/ira may stop about for times. In the 
case when a 0.02 p N cargo is applied the processivity re- 
duces to 64 steps. Finally, the load term fd/kBT can 
be viewed as the ratio between the work done by the 
kinesin to move such cargo and the thermal energy avail- 
able from surroundings (intracellular space), along the 
overall reaction. Since the load term is positive and the 
linear equilibrium term D eq is negative, we can conclude 
that the presence of a cargo reduces the processivity of 
motor by stopping the reaction before with respect to the 
case without load. 

Numerical values for the free-energy barriers were es- 
timated with Eq. (28) and are presented in Table IV, 
where we quantify the energy consumption of kinesin's 
translational motion for different times (steps). The 
change of sign of AG^ (at step 62) shows that the state 
MKT' becomes less favorable than MKT. This effect 
is also present with load (Fig. 8) at an earlier step. The 
increase of the energy barriers and G a -p is clear. 

VI. CONCLUSIONS 

In this work we proposed two biochemical models pro- 
viding the kinetic and energetic explanations of the pro- 
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Onm 



8nm 



16nm 



24nm 




Steps 



FIG. 8. Free energy landscape obtained via the Fourier ex- 
pansion Eq. (28) with k — 25 modes for a catalytic reaction 
of many cycles with cargo / = 0.02pN. The equilibrium val- 
ues of the mass fractions were n\ 
and n e ^ 

displacement, b) step 20 to 23. c) step 61 to 64. The values 
of the constants used are given in the caption of Figure 6. 



= 0.02, n T q = 0.78 
0.2. a) from cycle 1 to 3 and corresponding nm 



cessivity dynamics of kinesin, myosin and dinein-type 
molecular motors. Our approach is based on a well 
known model describing kinesin dynamics and consid- 
ers the presence of a competitive inhibition reaction by 
ADP. The two models discussed constitute a powerful 
tool to understand and describe in quantitative form the 
dynamics of translational molecular motors and therefore 
may be useful to give more precise descriptions of some 
cellular processes mediated by kinesins or myosins such 
as, for instance, exocytosis and endocytosis. 

We provided a new analytical procedure to reconstruct 
a continuous free-energy landscape of the cycle catalyst 
process starting from the corresponding biochemical re- 
action model. The obtained free-energy landscape is 
valid for a single motor and allows one to predict the 
heigh of the main free-energy barriers associated to the 
motion and the total number of steps given by the molec- 
ular motor for given physical and chemical conditions of 
the surroundings. That is the processivity of the motor 
is explained by means of energetic considerations. 

In addition, a collective description of the dynamics 
of translational molecular motors is also provided that 
allows to determine an analytical expression for the asso- 
ciated translational velocity as well as an expression for 
the stopping time of the molecular motors in terms of 



time and ATP concentration. 

Motor's processivity is discussed in quantitative form 
by using experimental data. The number of steps we pre- 
dict (60-66) agrees well with experimental observations 
and provides important quantitative information of the 
energy consumption during the process. In similar way, 
the average velocity and the time duration of a collective 
process is estimated via the kinetic description and also 
agrees (t s ~ 400s) with experimental reports on secre- 
tion processes mediated by kinesins and myosins. In this 
respect, the time dependence of motors' velocity may be 
very important since it may be used in coarse-grained 
models describing intracellular transport 20 . 
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Appendix A: Derivation of equations (9)-(ll) 

The main supposition in order to derive Eqs. (9) and 
(11), is to assume that the evolution of the coupled reac- 
tions (1) and (2) can be described by using a single degree 
of reaction £. This assumption is plausible because reac- 
tion (2) is in equilibrium, and therefore it allows to relate 
MKT with MKT' concentrations through the reaction 
constant . The main consequence of this assumption 
is that a fourth order polynomial in £ is recovered for 
the Gibbs free-energy associated to the enzyme reaction 
scheme. 

Let us start by expanding the free-energy change of 
the reaction (1) and using that elementary reactions obey 
the following relation between the chemical affinity Aj of 
j-th reaction and the corresponding degree of reaction: 
dGj = —Ajd^j. Then, it follows that for the reaction (1) 
the change on Gibbs free-energy is 



dGi = - (vmk^mKx + ^Ti^Ti H 



(Al) 



where the stoichiometric coefficients of reaction (1) are: 
vmK\ = 1, v Tl = 1 and vmkt x = — 1- Using now 
an ideal approximation for all chemical potentials (/i = 
ksThi \n\) we get 



dGi = -RT In 



vmk 1 
n MK 




vmkt x 
n MKT 


n MK eq 


n T eq 


n MKT eq 



d£i . (A2) 



Now, since we have assumed that reactions (1) and 
(2) are coupled, then it follows that at every time the 
concentration of the MKT complex can be calculated 

as: timkt = n° MKT + fi, or n M KT = n° MKT - f 2 , with 
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£2 the degree of reaction of the isomerization reaction. 
Comparing these two expressions we conclude that 



6 



(A3) 



The above relation is fundamental in the calculation of 
the free-energy Eq. (11) since it demonstrates that both 
reactions can occur with a same degree of reaction £. 

Now if we use equilibrium approximation for Eq. (2), 
we know that the equilibrium constant is defined as 29 



n 



n. 



Kl 



(A4) 



where V{ denotes the stoichiometric coefficient for reac- 
tants and v* the stoichiometric coefficient for products. 
For simplicity, in the following we will consider that all 
v* coefficients are negative. Hence, using this definition 
in reaction (2) we get 



K 



"mkt 2 _ n MKT 
MKT - K ] 



(A5) 



where the MifT-complex acts as the reactant and the 
complex MKT' as the product: vmkt 2 = 1 and 
v mkt! 1 — ~ 1- Noticing that for the reaction (1): 
vmkt x — — 1, whereas that for the reaction (2): vmkt 2 — 
1, then by using the result of (A3) we can establish the 
following relation 

(n°M KT + 6) I " M * T > 1 = (n° MKT - 6)"™ , (A6) 
and, with Eq. (A5) and (A6) we can also establish 

(n°MKT' + & f MKT i l = (n° MKT , - £1)—-' , (A7) 

where vmkt' = 1- Finally, from (A5)-(A7) it follows 
that 



VMKT X _( O 1 £■ \VMKT X 

MKT — \ n MKT "T Si; 



MKT/ 



■hY 



(A8) 

Eq. (A8) can be substituted into Eq. (9) in order to 
obtain Eq. (11), that is the effective Gibbs free-energy 
change for the coupled reactions (1) and (2) expressed 
in terms of the activated complex MKT 1 ', that is the 
measurable quantity 9 . 



Appendix B: Rescaling the initial mass fractions after each 
cycle 

During the advance of the reaction along the global 
reaction coordinate ADP is produced and ATP con- 
centration decreases by following the relation n^(^ g ) = 
Tiji — k cat T^ gi which follows directly from Eq. (22). Since 
there is a parallel inhibition reaction between the MK 
complex and ADP (ADP competes for the same active 
site of ATP), this reaction do not allows the formation of 



both MKT and MKT' complexes. The effect of inhibi- 
tion increases with £ g since the system is producing ADP. 
This decreasing behavior cannot be directly inferred from 
a kinetic equation, since the corresponding reaction ve- 
locity does not depends on the concentration of the com- 
plex MKT' . However we can make the assumption that 
the initial mass fraction n° MKT/ decreases in linear form 



with £ g , that is, n' 



MKT/\^9 



(&) 



n 



MKT/ 



K ig. This de- 
pendence is similar to that of n^(^), that follows from 
transition time theory 39 . Taking into account these con- 
siderations, we see that the Fourier coefficients depends 
on £ g in such a way that Eq. (11) transforms into 



dG 
RT 



{[n°MK ~ Z][n°AQ ~ £][n°MKT,(Z 9 ) ~ t] 



i eq n eq n eq 1 
MK T MKT/1 



f5/k B T}d£. 



(Bl) 



Starting from Eq. (Bl), we perform an integration over £ 
as in Section 3. The Fourier series can also be performed 
in which the Fourier coefficients are now given by the 
relations 

Aife) =/ cos(27r<)/ (n-x)(n-£g-x)(n-tg-x)dxd£, 
Jo Jo 

(B2) 



and 



B n (Q =/sin(27mO 



(n- 



X)(n-ig-x)(n-ig-x)dxdi. 

(B3) 

From Eqs. (B2) and (B3) it follows that, after performing 
integrations, linear and quadratic terms in ^ g will remain. 
Thus, A n {i g ) and B n (^ g ) will be second order polynomi- 
als in and not constants as in Eq. (14). For notation's 
simplicity, we dropped out explicit dependencies on k cat 
and K Q in the last equations. 



Appendix C: Effects of non competitive inhibition by ADP 
and inhibition by orthophosphate. 

As mentioned in the text, one might also consider non- 
competitive inhibition by ADP, competitive and non- 
competitive inhibitions by orthophosphate and also a 
mix ADPoP^ non-competitive and competitive inhibi- 
tions. These effects will introduce new inhibition con- 
stants: Kjj for non-competitive ADP inhibition, K l 

p 

for competitive Pi inhibition, K IT for non-competitive 

7-> ADPoP ADPoP . 

Pi inhibition and K IT and K T lor the mixed 
cases. Considering these effects the numerator and de- 
nominator of Eq. (20) become 



dk c 



and 



1 



[D] , [P] , [D][P] 



K 



11 



K 



11 



K 



11 



[D][P] 

j y ADPoP 

K I 



(CI) 



(C2) 
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respectively. The above expressions together with the 
time dependence of the concentration of orthophosphate 
[P]t leads to a new expression for ^ a£C and h t in Eq. (24). 
However, as it was previously indicated, the effect of 
these type of inhibitions is significative only at low ATP 
concentrations, that is, only in a time interval near the 
end of the process 9 . 
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